Quantum State Reconstruction of Many Body System 
Based on Complete Set of Quantum Correlations Reduced by Symmetry 
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We propose and study a universal approach for the reconstruction of quantum states of many body 
systems from symmetry analysis. The concept of minimal complete set of quantum correlation 
functions (MCSQCF) is introduced to describe the state reconstruction. As an experimentally 
feasible physical object, the MCSQCF is mathematically defined through the minimal complete 
subspace of observables determined by the symmetry of quantum states under consideration. An 
example with broken symmetry is analyzed in detail to illustrate the idea. 
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Introduction: The concept of quantum state is funda- 
mental in quantum physics. Quantum state can provide 
us with a complete knowledge of the considered system 
to predict the possible results of any measurement on the 
system m . But the puzzling property of quantum coher- 
ence, which is intrinsic to quantum state, imposes insur- 
mountable limitations on our ability to fully reconstruct 
quantum state of a single quantum system by devising a 
complete set of measurements on the system 2j. 

Actually, we cannot make a measurement on a sin- 
gle quantum system without back-action and the non- 
cloning theorem forbids us from producing an exact 
copy of some quantum system in an unknown state. But 
on the other hand, if identical copies of quantum sys- 
tem, or an ensemble, is initially prepared in the same 
quantum state, it is possible to estimate the unknown 
state of the quantum system by carrying out appropriate 
measurements on each copy 4]. This is just the idea of 
quantum tomography, which has recently become a fash- 
ion in quantum information physics, motivating many 
theoretical and experimental investigations . 

We notice that most works on quantum tomography 
or other quantum state reconstruction methods focus 
on continuous variable system and single qubit system. 
In this letter we will deal with the many particle case, 
which possesses distinguished features such as quantum 
entanglements . It has become well known that quan- 
tum correlation and the corresponding Green function 
are important conceptual tools to probe the non-local 
nature of many body system. On the other hand, in 
the study of quantum information, quantum entangle- 
ment turns out to be a crucial concept and thus becomes 
a fundamental element in the research about the non- 
locality of quantum system. Hence, it is doubtlessly im- 
portant and meaningful to probe the possible connection 
between quantum entanglement and quantum correlation 
0. Along this line, some progress has been made. For 
example, for the spin chain system with the symmetry of 
SO{2), a connection has been revealed between the con- 
currence, which is a measure of quantum entanglement. 



and the correlation functions of first order and second 
order, which is supposed to reflect the property of cor- 
relation |8j. This connection suggests a possible connec- 
tion between quantum entanglement and quantum phase 
transition and prompts us to introduce in this 

letter a general framework for studying such connection. 

Our starting point is the basic observation that com- 
plete information about quantum state of many body 
system can be obtained from the expectation values of 
properly chosen single body observables and their corre- 
lation functions of various orders ^'^^ our idea is 
to reduce the number of times of measurement needed to 
determine a quantum state by making full use of the sym- 
metry that the state is supposed to possess. The main 
purpose of this letter is to make clear the mechanism of 
this reduction and present interesting examples to illus- 
trate it. To this end we will introduce the minimal com- 
plete subspace of observables (MCSO) associated with a 
symmetry group, which mathematically determines the 
quantum state with this symmetry, and the minimal com- 
plete set of quantum correlation functions (MCSQCF) 
associated with the MCSO, which determines the quan- 
tum state physically in some sense. 

State Reconstruction Based on Correlation Functions: 
Throughout this letter, let V = =Wi®W2® 
■■■(g) Wm be the Hilbert space of a many spin system. 
Here the Wi 's stand for the state space W of a, spin 
system. For an arbitrary finite dimensional Hilbert space 
S we denote by A{S) the set of Hermitian operators on S. 
It is well known that A{S) is an Euclidean space with the 
inner product ( , ) defined as {A, B) = iv{AB), ^ A, B e 
A{S), where tr is the trace over W. 

Now let p £ A{V) be a density operator and {Ai} a 
basis of A{V) . Then p can be written as a linear combina- 
tion of Ai's: p — J2i CiAi. We notice that the coefficients 
Ci's are determined by the expectation values {tr(pAi)} 
of the observable At in the state p. In particular, if {At] 
is an orthonormal basis then we have Ci = ivpAi. Obvi- 
ously, from mathematical point of view it is trivially true 
that p is completely determined by {{p,Ai)}. But phys- 
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ically it means the possibility of determining a quantum 
state by measuring some properly chosen observables. 

To be precise, let us choose a basis {Oi} of A{W). The 
Oi's can be understood as observables of single particle 
system. Then {(^^^lOi^} is a basis of A{V), where each 
Oi- belongs to {Oi}. Hence we can identify the state p 
by the quantities 

tr(p 0,J ^ (®r=iO,Jp, (1) 

which are just correlation functions of the considered sys- 
tem. 

For spin 1/2 system interesting things happenj3,l3- In 
this case we can take {Oi} to be {cto: fi, cr2, tra}, where 
(To is the identity operator and cri,(T2,(73 are the Pauli 
operators Ux , Cj, , CTz . Let us introduce the standard corre- 
lation functions Gf^i^.-.i^ — (®fe=iO'ifc)p, associated with 
the state p. Clearly, from both theoretical and experi- 
mental points of view standard correlation functions are 
accessible physical quantities. According to the discus- 
sion above, for spin 1/2 system a state is determined by 
these ideal quantities. 

We are not so lucky when we consider higher spin sys- 
tem. The main reason is, when the dimension of W is 
larger than 2 it is impossible to span ^(W^) with just 
the spin operators and identity operators. This forces us 
to consider the correlation functions not related to the 
spin operators in a direct way. From experimental point 
of view such correlation functions are less physical than 
the standard correlation functions. But from theoreti- 
cal point of view there is no essential difference between 
them. 

Now a natural question is: how many correlation func- 
tions are needed to determine a quantum state? Ob- 
viously, if we have no knowledge about the state, then 
the number of correlation functions needed is exactly the 
dimension of A{V). But if we are sure that the state pos- 
sesses some symmetry, then it turns out that the number 
might be greatly reduced. Let us elucidate this point in 
detail as follows. 

Reduction to the Minimal Complete Subspace of Ob- 
servables by Symmetry: Now, we discuss symmetry of 
density operator and the role symmetry can play in the 
reconstruction of density operator. 

The symmetry group of the density operator p e A{V) 
is defined to be the subgroup Up of U{n) with the prop- 
erty upu'' = p, y u G Up. If [/ is a subgroup of Up we 
will say that p possesses the symmetry of U. In practice, 
it is highly possible that only partial information about 
the symmetry of a density operator is available. In other 
words, one may only be sure that some group is a sub- 
group of the symmetry group of a density operator. This 
being the case, it is desirable to know how to identify 
a density operator among a family of density operators 
whose symmetry groups contain a specific unitary group 
U (a subgroup of U{n) by definition). This is the idea 
underlying the following discussion. 



A symmetry of U determines a family Ou of observ- 
ables in the following way: A € Ou if and only if 
uAu^ ^ A, y u e U. By definition, if A £ Ou then 
its symmetry group contains U as a subgroup and vice 
versa. Thus p belongs to Ou if and only if it possesses 
the symmetry of U. Clearly, Ou is a subspace of A{V). 
We call it the MCSO associated with U. 

As A{V) is an inner product space it can be decom- 
posed as A{V) — Ou®Oij, where O^ is the orthonormal 
complement to Ou- Let {Ai \ be an arbitrary basis of Ou- 
Then each element A of Ou is uniquely determined by 
the quantities {{A^Ai)}. This justifies calling Ou com- 
plete. Consequently, if p possesses the symmetry of U 
then the number of measurements needed to determine 
p is just the dimension of the MCSO associated with U- 

Ncxt let us determine the structure of Ou- 

Proposition 1. If U is connected, then iOu is the 
centralizer of the Lie algebra of U in u{n), the Lie al- 
gebra of U{n), where n is the dimension of V- ( Here 

Proof. We notice that iA{V) is exactly the Lie alge- 
bra of U{n). Let i A be an arbitrary element of the Lie 
algebra of U and i B an arbitrary element of i Ou- Then 
we have iB = exp{it A){i B) exp{—it A) or 

iB:^{Ade''^){iB)=e'''^'^{iB). (2) 

It follows that iOu is a subset of the centralizer of the Lie 
algebra of U- That they are actually identical is now a 
consequence of the fact that if U is connected, then it can 
be generated by 1-parameter subgroups. This completes 
the proof of the proposition. 

When [/ is a 1-parameter subgroup the structure of 
Ou can be given a concise description. Let us consider 
this case in more detail. Suppose that U = exp {—it A) 
where A G A{V), then according to Proposition 1, an el- 
ement B belongs to Ou if and only if [A, B] — 0. As A is 
an Hermitian operator there exists an orthonormal basis 
7 of 1^ with respect to which it has the matrix represen- 
tation of block diagonal form (the meaning of notation 
here should be self-evident): A = diag{v4i , , • • • , A„}, 
where Ai = Xilm, K being an eigenvalue of A and /„; 
is the unit matrix of rank n^, being the multiplicity 
of A^. From this observation, the following proposition 
follows directly. 

Proposition 2. If U is connected then an element B 
lies in Ou if and only if its matrix representation takes 
the form B — diag{Bi, B2, ■ ■ ■ ,-Bn}, with respect to the 
same basis 7 defined above, where Bi is an Ui by Ui 
matrix. 

Symmetry of Reduced Density Operators: In this sec- 
tion we consider symmetry of the so called reduced den- 
sity operator, which is derived from the density operator 
after a process of taking partial trace. We are interested 
in the symmetry that survives such process. 

Let V = Vi ®V2 and p a density operator in A{V). 
The reduced density operator tr2p is defined to be an 
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element of ^(Vi) satisfying the relations (xi, {tr2p)x2) = 
J2i{^i^yi^ P^2^yi)y G Vi, where {j/j} is a basis of 

V2. Note that we have used the same notation ( , ) for the 
inner products in the different spaces. It is easily check 
that tr2p is also a density operator. Now the following 
lemma comes as a direct consequence of the definitions. 
We would rather omit the proof. 

Lemma 1. If p possesses the symmetry of Ui®U2 
then tr2p possesses the symmetry of U\. 

Although the idea of studying how symmetry survives 
the reduction process might seem trivial, there are non- 
trivial examples to support it. Let us consider two models 
of such examples. For convenience, from now on we will 
use the Dirac notation. 

Suppose that A e A{V) has the form A = Ai 1 + 
1(^)^2, where Ai e A{Vi) and A2 € A{V2). Let {ip) be a 
normalized eigenvector of A. Then it is easily seen that 
the density operator p = \ip) (ipl possesses the symme- 
try of the 1-parameter subgroup cxp [it A) that can be 
written as exp {it A) = exp {it Ai) (g) exp {it A2). Thus the 
reduced density operator tr2/9 possesses the symmetry of 
exp {it Ai) according to Lemma 1. 

In the special case that p has the symmetry of 1- 
parameter subgroup there is another mechanism of pass- 
ing symmetry to the reduced density operator, which is 
different from the mechanism shown above. Let A be 
an Hermitian operator of the form A = Ai ^ A2 where 
Ai € A(Vi) and A2 G A{V2) and a normalized eigen- 
vector of A. Then for the density operator p = we 
have the following result. 

Lemma 2. The reduced density operator tr2P pos- 
sesses the symmetry of exp(ityli) if the eigenvalue cor- 
responding to lip) is nonzero. 

Proof. Let be an orthonormal basis of Vi con- 

sisting of the eigenvectors of Ai and an orthonor- 

mal basis of V2 consisting of the eigenvectors of A2. Sup- 
pose that Ai\(f>i) = Xi\(f>i) and A2\Lpj) = fij\Lpj). Since 
lip) is an eigenvector of A there exist and pi such that 
A\ijj) = XkPi\ijj). By definition, for arbitrary two elements 
\(t>i-^) and {(pi^) of the basis {If/ii)} wc have {(pi-^^ \tT2p\(t>i2) = 

,(pj), where | (pi^ , 
It follows that {(f)i^\tr2p\4'i2) = Oi unless Aj^ = Aj^. On 
the other hand, we have 

J' = e'*(^--^-)(</>iJtr2p|<^i,). (3) 

for P = {<l)i^\exp{it Ai)tT2pexp{— it Ai)\(j)i2). This just 
means that 

p = I i^ii I tl'2P 10*2), if Ail = /^x 

1^ 0, otherwise. ^ ' 

Thus we conclude that exp (it Ai) tr2p exp (— if ^i) = 
tr2/3 and complete the proof. 

Minimal Complete Set of Correlation Functions: In 
this section the problem of reconstruction of quantum 
state is reduced to the problem of determination of a set 
of correlation functions. 



Let {O,} be a basis of A{W). Then {«)^^iO,J is a 
basis of A{V). We call it the measurement basis of A{V) 
associated with {Oi}. A subset of the measurement ba- 
sis of A{V) is called a reconstruction basis of the MCSO 
associated with {Oi} if each element of the MCSO can 
be written as a linear combination of elements in the 
subset and a reconstruction basis is called minimal if it 
does not contain any smaller reconstruction basis. Note 
that according to this definition elements of a reconstruc- 
tion basis of the MCSO do not necessarily belong to the 
MCSO. 

Proposition 3. There exists a unique minimal recon- 
struction basis of an MCSO associated with each {Oi}. 

Proof. Choose an arbitrary basis (3 of the considered 
MCSO and define a subset 71 of the measurement basis 
{(g^iOifc} as follows: an element lies in 71 if and only if 
it appears as a non zero term in the expression of some 
element of f3 in terms of {^feLiOifc }. Clearly 71 is a 
minimal reconstruction basis associated with {Oi}. We 
claim that it is in fact the unique minimal reconstruction 
basis of the considered MCSO associated with {O,}. 

Let 72 be another minimal reconstruction basis asso- 
ciated with {Oi}. Then each element of (3 can be ex- 
pressed as a linear combination of the elements of 72 or 
as a linear combination of the elements of 71. Thus by 
equating the two expressions we conclude that 71 C 72, 
considering the fact that both 71 and 72 are subsets of 
the measurement basis {^/JLiOi^}. That 71 and 72 are 
actually identical now follows from the minimality of 72. 
This completes the proof of the proposition. 

We are now in a position to introduce the concept 

of MCSQCF. Let Ou be an MCSO and {®r=i0jj 
the minimal reconstruction basis of Ojj associated with 
{Oj}. Please note the notational difference between 
{®feLiOjfc}and {(XifeLiOifc}. Here the former is used to 
denote a subset of the latter. With respect to {(g^^Oij, }, 
the MCSQCF associated with Ou is defined to be the set 
{{'^T=i^jk) p} °f quantum correlation functions where p 
stands for a density operator in Ou- The next proposi- 
tion shows that each p in Ou is completely determined by 
{{'^T=i^jk) p} b^it not by any proper subset of it. This 
is just what "complete" and "minimal" mean in the def- 
inition of MCSQCF. 

Proposition 4. A density operator p in Ou can be 

expressed in terms of the MCSQCF {(®5!LiOj Jp}, but 
cannot be expressed in terms of any proper subset of 
{(0^iOjj.)p}, regardless of the basis of V on which p 
acts. 

Proof. Let {Bi] be an orthonormal basis of Ou- We 
then have p = Y^-{p,Bi)Bi = ^.(trpi?i)i?i, where each 
coefficient trpSi can be written as a linear combination of 
{{^'k=iOjk)p}- Moreover, each ((X)^^iOj^)p will appear 
in some ivpBi by the minimality of {®k=i^3k}- Hence p 
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can be written as 

P= (®r=iOjJpi?j,j,...,„, (5) 

where Bj-^j^...j^ is a linear combination of {Bi}, which 
cannot be a zero matrix with respect to any basis of 
V since {Bi} is hnearly independent. This proves the 
proposition. 

Example with Broken Symmetry : Now let us present a 
concrete example to illustrate the results obtained above. 
We consider the case that W is the state space of spin 1/2 
particle. Let | t) and | J,) be the normalized eigenvectors 
of the Pauli operators CTz and take V = Vi (E) V2,Vi = 

Wiig)W2,V2 = W3(g)---<E Wm- 

We suppose the state p to be measured possess the 
symmetry generated hy A — (X'^'Lif^, namely the sym- 
metry of exp (it A) , where m is even. Then according 
to Lemma 2, the reduced density operator p = tr2p 
possesses the symmetry of J7 = exp{itAi) with Ai = 
az® (Jz- Consequently, by Proposition 2, with respect to 
the basis {| tt), I Ti),| it)}, the reduced density 

operator is of the 8- vertex form 

~P = , (6) 

\ ri ) 

where the r^'s are real numbers with the restriction 
— 1 and z* stands for the complex conjugation 
of Zi. One can check that {I®!, 1 ® az, <7z ax ® 
o'x, o'y (8> ay, az az, a^ ® ay, ay ® a^) is the unique 
reconstruction basis of Ou associated with {ai\ and that 
{Gqq, Gpg, GgQ, C-^^, G22, G33, G^2, G21} is the cprre- 
sponding MCSQCF. The explicit expression of p in terms 
of the MCSQCF is as follows: ri = (G^g + G|o + G^g + 
l)/4, r2 = (Gg3 - G§o - Gg3 + l)/4, r, = (Gg„ - Gi3 - 
Ggg + l)/4, r4 = (Gg3 - Gg3 - Ggo + l)/4, z, - (G?i - 
G^2 + *G?2 + *G^i)/4, Z2 = (G?i + G^2 + zGf^, - iG^^)!^. 

This example has a non-trivial connection with the 
Ising model with the Hamiltonian 

m — 1 m 

H = -jY,aiai^^-rgY,ai, J>^. (7) 

When the external field g = 0, this model possesses the 
S'0(2) ® Z2 symmetry where Z2 represents the rotations 
of all spins by 180 degrees about the x axis while 50(2) 
represents the rotations of all spins by arbitrary degrees 
about the z axis. However, this system has two degener- 
ate ferromagnetically ordered ground states | f)®™ and 
I I)®™, which break the discrete symmetry Z2. Obvi- 
ously, as a superposition of these two ground states, the 
cat statelV') = a| T)®" + fi\ i)®" possess the symmetry 
of exp [it A) mentioned above since it is the eigenvector 
oi A = (Eijli^z corresponding to the eigenvalue 1, and 



then the unknown a and /3 satisfying |ap -|- = 1 can 
be determined by the above quantum state reconstruc- 
tion method. When g 7^ 0, the reduced density operator 
p = tr2p derived from the ground state or the thermal 
equilibrium state(a function of hamiltonian) is of the 6- 
vertex form for the existence of the total spin projection 
a],. This result has appeared in many places, e.g., in 
the references 

Concluding remarks : Based on the viewpoint of sym- 
metry, we have proposed a quantum state reconstruc- 
tion method for density operators or reduced density op- 
erators of the many body system. Our central idea is 
to utilize the MCSQCF determined by the symmetry of 
quantum states rather than the symmetry of the Hamil- 
tonian. The present work might be regarded as a general- 
ization of the "particle creation-annihilation expression" 
for the reduced density matrices of identical particles in 
Ref.|l2l|. which reveals the essence of quantum conden- 
sation. Thus, we hope that the results in this letter can 
be directly applied to explore the quantum and classi- 
cal critical phenomenon for various quantum many body 
systems. 
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